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Abstract 

Using in-medium hadron properties according to the Brown- Rho scaling due to the chiral 
symmetry restoration at high densities and considering naturalness of the coupling constants, 
we have newly constructed several relativistic mean-field Lagrangians with chiral limits. The 
model parameters are adjusted such that the symmetric part of the resulting equation of state 
at supra-normal densities is consistent with that required by the collective flow data from high 
energy heavy-ion reactions, while the resulting density dependence of the symmetry energy 
at sub-saturation densities agrees with that extracted from the recent isospin diffusion data 
from intermediate energy heavy-ion reactions. The resulting equations of state have the 
special feature of being soft at intermediate densities but stiff at high densities naturally. 
With these constrained equations of state, it is found that the radius of a 1.4M Q canonical 
neutron star is in the range of 11.9 km<R<13.1 km, and the maximum neutron star mass 
is around 2.OM close to the recent observations. 

Keywords: Equation of state, nuclear matter, symmetry energy, relativistic mean-field 
models, chiral limits. PACS numbers: 21.65. +f, 26.60.+C, 11.30.Rd 

1 Introduction 

The Equation of State (EOS) of isospin asymmetric nuclear matter plays a crucial role in many 
important issues in astrophysics, see, e.g., Refs. [1, 2, 3]. It is also important for understanding 
both the structure of exotic nuclei and the reaction dynamics of heavy-ion collisions, see, e.g., 
Ref. [4]. Within the parabolic approximation, the energy per nucleon in isospin asymmetric 
nuclear matter can be written as E/A = e(p) + E sym (p)5 2 where e(p) is the EOS of symmetric 
nuclear matter, the E sym (p) is the symmetry energy and 5 = (p n —Pp) / P is the isospin asymmetry. 
Both the e(p) and the E sym (p) are important in astrophysics although maybe for different issues. 
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For instance, the maximum mass of neutron stars is mainly determined by the EOS of symmetric 
nuclear matter e{p) while the radii and cooling mechanisms of neutron stars are determined 
instead mainly by the symmetry energy E sym [l, 5]. The nuclear physics community has been 
trying to constrain the EOS of symmetric nuclear matter using terrestrial nuclear experiments 
for more than three decades, see, e.g. [6], for a review. On the other hand, a similarly systematic 
and sophisticated study on the density dependence of the symmetry energy E sym using heavy-ion 
reactions only started about ten years ago stimulated mostly by the progress and availability of 
radioactive beams [7]. Compared to our current knowledge about the EOS of symmetric nuclear 
matter, the symmetry energy E sym is still poorly known especially at supra-normal densities [8] 
given the recent progress in constraining it at densities less than about 1.2po using the isospin 
diffusion data from heavy-ion reactions [9, 10, 11]. 

The aim of this work is to investigate the EOS of isospin asymmetric nuclear matter within 
the relativistic mean-field (RMF) model with in-medium hadron properties governed by the BR 
scaling. From the point of view of hadronic field theories, the symmetry energy is governed 
by the isovector meson exchange. Studying in-medium properties of isovector mesons is thus 
of critical importance for understanding the density dependence of the symmetry energy. We 
first construct model Lagrangians respecting the chiral symmetry restoration at high densities. 
The model parameters are adjusted such that the symmetric part of the resulting EOS at 
supra-normal densities is consistent with that required by the collective flow data from high 
energy heavy-ion reactions [6], while the resulting density dependence of the symmetry energy 
at sub-saturation densities agrees with that extracted from the recent isospin diffusion data 
from intermediate energy heavy-ion reactions [9, 10, 11]. The constrained EOS is then used to 
investigate several global properties of neutron stars. 

2 Relativistic mean-field models with chiral limits 

In-medium properties of the isovector meson p can be studied through the special symmetry 
breaking and restoration. The local isospin symmetry in the Yang-Mills field theory, where the 
p meson may be introduced as a gauge boson of the strong interaction, can serve as a possible 
candidate to study the in-medium properties of p meson. However, since ttN interactions actually 
dominate the strong interaction in hadron phase, it was rather difficult to understand how the 
in-medium properties of the massive p meson could be consistent with the restoration of local 



2 



isospin symmetry. On the other hand, within the microscopic theory for the strong interaction, 
namely, the QCD which is a color SU(3) gauge theory, the chiral symmetry is approximately 
conserved. The spontaneous chiral symmetry breaking and its restoration can be manifested 
in effective QCD models. Based on the latter, Brown and Rho (BR) proposed the in-medium 
scaling law [12] implying that hadron masses and meson coupling constants in the Welacka 
model [13] approach zero at the chiral limit. The scaling was treated in the hadronic phase 
before the chiral symmetry restoration. 

As an effective QCD field theory, the hidden local symmetry theory has been developed 
to include the p meson in addition to the pion in the framework of the chiral perturbative 
theory by Harada and Yamawaki [14, 15] and it is shown that the p meson becomes massless 
at the chiral limit. This supports the mass dropping senario of the BR scaling. There are 
also experimental indication for the mass dropping, i.e., the dielectron mass spectra observed 
at the CERN SPS [16, 17], the to meson mass shift measured at the KEK [18] and the ELSA- 
Bonn [19], as well as the analysis of the STAR data [20, 21]. However, data from the NA60 
Collaboration for the dimuon spectrum [22] seem to favor the explanation of p meson broadening 
based on a many-body approach [23]. So far, the controversy is still unsettled [21]. The chiral 
symmetry and its spontaneous breaking are closely related to the mass acquisition and dropping 
of hadrons. Since the chiral symmetry is a characteristic of the strong interaction within the 
QCD, it is favorable to include in the RMF models effects of the chiral symmetry through 
the BR scaling law. However, this does not mean that the contribution of the many-body 
correlations [24] is excluded. Actually, the contribution of the many-body correlations can be 
included phenomenologically into the RMF models to reproduce the saturation properties of 
nuclear matter. 

The in-medium p meson plays an important role in modifying the density dependence of the 
symmetry energy. For most RMF models that the p meson mass is not modified by the medium, 
the symmetry energy is almost linear in density. The introduction of the isoscalar-isovector 
coupling in RMF models can soften the symmetry energy at high densities [3]. Meanwhile, 
it reproduces the neutron-skin thickness in 208 Pb as that given by the non-relativistic models 
(about 0.22fm) [25], consistent with the available data [26]. In well-fitted RMF models that 
give a value of incompressibility k = 230 MeV, a large coefficient of non-linear self-interacting 
u) meson term is required [25] and thus the naturalness breaks down. Moreover, the isoscalar- 
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isovector coupling in the RMF models increases the effective p meson mass with density, which 
leads the model to be far away from the chiral limit. 

The Walecka model with the density-dependent parameters is the simplest version to incor- 
porate the effects of chiral symmetry. The Lagrangian is written as 

C = ^[i 1 ^-M*+g:a-g:^-g;^T 3 b^ + ^(d,ad^a-m: 2 a 2 ) 

-\f»„F^ + \m*J^ - \B, V B^ + \mfb,X U) 

where ip,a,uj, and 60 are the fields of the nucleon, scalar, vector, and isovector-vector mesons, 
with their in- medium scaled masses M*,m*,m*, and m*, respectively. The meson coupling 
constants and hadron masses with asterisks denote the density dependence, given by the BR 
scaling. The energy density and pressure read, respectively, 

2 f k Fi 



£ = \CIP 2 + \cyS 2 + \cl{m* N - M*) 2 + E 7A3 / ^ (2) 

i=p,n v / 

P = + \c 2 p pH 2 - \cl{m* N - M*) 2 - Sop + \ £ /*"cft ^ (3) 



i=p,n 



where C w = gZ/m*, C p = g*/m* p , C a = m* a /g*, E* = y'k 2 + m* N 2 with m* N = M* — g*a the 

effective mass of nucleon, and kp is the Fermi momentum. The incompressibility of symmetric 

matter can be expressed explicitly as [27] 

^ dp n ^ dC w dE F <9£ 

K = gpi t = M Cl + 2C^— + — -—) (4) 

where the chemical potential is given by p = d£ /dp and the Fermi energy is Ep = \Jkp + m* N 2 . 

The rearrangement term is essential for the thermodynamic consistency to derive the pressure 

in (3) and its expectation value in the mean field So is given by 

2 dC u 2 , 2r dC p ~8C , 2 dM* 

So = ~P C„— - p 5 C p — - C„ — (m N - M ) - p s — . (5) 

The density dependence of parameters is usually described by the scaling functions that are the 
ratios of the in-medium parameters to those in the free space. The choice of scaling functions and 
their coefficients are constrained by the saturation properties of nuclear matter and experimental 
data about the in-medium mass dropping of vector mesons. Moreover, we also use as a constraint 
the pressure within the density range 2 — 4.Gpo extracted from measurements of nuclear collective 
flows in heavy-ion collisions [6]. The scaling function may take the form [28]: 



with y = 0.28 for the vector meson mass, giving $>(po) = 0.78 found in QCD sum rules [29]. 
Recently, in Ref. [30] where the memory effect in dimuon yield was studied by considering the 
mass dropping of p meson, the authors cited a scaling function [31] 

Hp) = i - yp/po (7) 

with y = 0.15. Data extracted from the j — A reaction by the TAP collaboration indicate a value 
of y ps 0.13 [32]. In addition, data from the KEK photon-induced nuclear reaction indicate that 
the u> meson mass dropping is about the same order of magnitude at the saturation density [33]. 

Song firstly built effective models based on the BR scaling using the scaling functions (6) 
for both hadron masses and vector coupling constants [28]. A reasonable incompressibility for 
nuclear matter was obtained by introducing the non-linear self-interacting meson terms with 
coefficients satisfying the hypothesis of naturalness that is originated from the chiral symmetry 
and QCD scaling [34, 35, 36, 37]. In Ref. [38], along the line of [28], the BR scaling function 
(6) was considered only for hadron masses with a small value of y in the RMF models to study 
nuclear matter properties. In a more recent work [39], the BR scaling function (6) was taken 
for the scalar and vector meson coupling constants with respective values of y, and the scaling 
function (7) was taken for hadron masses in the RMF models without the self-interacting meson 
interactions. Though the models built in these works can give rather good descriptions of nuclear 
saturation properties, the pressures calculated in the density region of p = 2 — 4.6po> however, 
are still far away from that extracted from measurements of nuclear collective flows in heavy-ion 
collisions [6]. People have already tried to improve the situation by including the non- linear 
meson self-interacting terms. Unfortunately, it is not satisfactory because even unreasonably 
large coefficients of non-linear meson terms that break down the hypothesis of naturalness can 
not reduce the pressure lower enough to pass the experimental pressure-density region. 

The effective nucleon mass is not dominated by the BR scaling (at least at the normal 
nuclear matter density) since it is usually around 0.65M* at normal density while the mass 
dropping given by the BR scaling is less than 15% (y < 0.15 in (7)). This implies that the 
scalar meson coupling constant that plays a crucial role in the effective nucleon mass can adopt 
a different density-dependent scaling from that of the vector meson. In particular, the high 
pressure predicted by various models at high densities should be lowered, and this requires the 
scalar coupling constant to decrease more slowly. Furthermore, if the scaling function (7) for 
the to meson mass is preferred by experiments, one may take the same scaling function (7) for 
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Table 1: Parameter sets fitted at the saturation density po = 0.16fm~ 3 . The vacuum hadron 
masses are M = 938MeV, m a = 500MeV, m w = 783MeV and m p = 770MeV except for m a = 
600MeV for the parameter set S3. The coupling constants given here are those at zero density. 
For parameter sets SL3 and S3, the non-linear a self-interacting coefficients are introduced (see 
text). The parameter set SL1* has two more parameters y p = 0.654 and y^ = 0.0365. The 
symmetry energy is fitted to 31.6MeV at p = 0.16fm~ 3 for all models. The critical density p c 
for the chiral symmetry restoration is given by the value y in (7) for zero hadron mass. 
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SL1 


8.6388 


10.4634 


3.7875 


0.126 


0.234 


230.0 


1.0 


0.679 


7.94 


SL1* 


9.7414 


12.5535 


5.8644 


0.126 


0.238 


230.0 


1.0 


0.600 


7.94 


SL2 


6.1664 


10.9682 


3.9866 


0.11 


0.381 


219.5 


0.89 


0.763 


9.09 


SL3 


9.8627 


12.4928 


3.6128 


0.126 




250.0 


1.0 


0.620 


7.94 


S3[28] 


5.3210 


15.3134 


3.6035 


0.28 




250.0 


0.78 


0.617 





the coupling constant of to meson to avoid the infinity of pressure at the chiral limit where the 
scalar density vanishes. In this way, the effect of density dependence from the vector meson part 
is cancelled out since the energy density (3) only relies on the ratios C u and C p . Generally, we 
may take the scaling functions for the coupling constants of vector mesons as 

^ , \ i - yp/po i x i - yp/po ( ^ 
i + y P p/po i + y^p/po 

For hadron masses (including nucleons, if they have), (7) is taken with the same value of y used 
in (8). For the a meson coupling constant, the same form as (6) is taken but with a coefficient 
denoted by x: 

3 Results and discussions 

We first adjust the parameters to reproduce the saturation properties including the binding 
energy per nucleon E/A — M = —16 MeV, the zero pressure, the incompressibility k and the 
effective nucleon mass m* N at saturation density po = 0.16fm~ 3 . The resulting parameter sets 
SL1 and SL2 and the corresponding saturation properties are tabulated in Table 1. In SL1 the 
nucleon mass scaling is not considered and the effective nucleon mass is just m* N = M — g*a. In 
SL2, the nucleon mass scaling is included, and a much larger effective nucleon mass at normal 
density is obtained. There are just two coefficients x and y used in the scaling functions in SL1 
and SL2. Without the inclusion of the non-linear meson self-interacting terms, the saturation 
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Figure 1: The pressure as a function of density for different models. The shaded region is given 
by experimental error bars [6]. 

properties at po = 0.16fm~ 3 (see Table 1) and the pressure within the density region p = 2— 4.6po 
are both nicely reproduced with the SL1 and SL2 (see Fig. 1). It is worth mentioning that the 
vector potential which is quadratic in density for the constant is known to result in a higher 
pressure above the experimental region. The softening of the pressure here is attributed to the 
contribution of rearrangement terms. For the SL1, only the rearrangement term from the a 
meson survives. 

The parameter set SL3 does not consider the scaling of scalar meson coupling constant but 
includes in the Lagrangian the non-linear a self-interacting terms U(a) = <?2C 3 /3 + g^a^/A with 
coefficients gi = 23.095 and g% = —29.678. Though these large coefficients are needed to fit the 
saturation properties, they seem to be inconsistent with the hypothesis of naturalness [35, 36]. 
In S3, we introduce a g<i = —1.096 to obtain the given k in Table 1. This is a little different 
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from the original one [28]. As a comparison, we can see in Fig.l that the SL3 and S3 parameter 
sets are not consistent with the pressure constrained by the collective flow data. 
With (3), the symmetry energy in the RMF models can be derived as 

_l d\£/p) _l 2 kl 

Esym ~ 2^2- ~ lf P P + g^- (10) 

It consists of contribution from the p meson (potential part) and nucleons (kinetic part). In 
SL1 and SL2, we adopt the same scaling function for coupling constants of p and u> mesons: 
^p(p) = ^l<j(p) = 1 — yp/po- in this way, the ratio C p is just a constant that does not rely on 
the density. An almost linear dependence of the symmetry energy on the density is expected 
from Eq.(10). The symmetry energy as a function of density is shown in Fig.2. In Refs. [10, 11], 
the symmetry energy extracted from the isospin diffusion data is parameterized as E sym = 
31.6(p/po) 7 w hh 0.69 < 7 < 1.05 . All the parameter sets based on the BR scaling discussed 
above lead to the symmetry energies between those parameterized with 7 = 0.69 and 1.05 in the 
whole density region. In Ref. [11], the authors pointed out that with the in- medium nucleon- 
nucleon cross sections, a symmetry energy of E sym (p) = 31. 6(/?/ po) ' 69 for p < 1.2po was found 
most acceptable compared to the isospin diffusion data. However, the symmetry energy at higher 
densities is not constrained at all. It is thus interesting to examine predictions within the RMF 
models with the chiral limit at high densities. 

Our strategy is to first fit the symmetry energies constrained at low densities by modifying 
the coefficient y p in (8), and then predict the symmetry energy at high densities. The modified 
symmetry energy is shown in Fig.2 with the dashed blue curve denoted as SL1*. Considering 
the nucleon effective mass at normal density is a little large in SL1, we introduce a coefficient 
y u in (8) to reduce it to the value of 0.6M. The new parameter set is called SL1* and also listed 
in Table 1. The different parameterizations in SL1 and SL1* lead to significantly difference in 
pressure at high densities as shown in Fig.l. However, the suppression of the symmetry energy 
in SL1* at high densities is dominated by the density dependence of the ratio C p in (10). This 
can be seen by comparing results shown in Table 1 and Fig.2 that the density dependence of 
the symmetry energy does not change much by the different nucleon effective masses of various 
parameter sets. In the density-dependent RMF model [24], the dropping of ratios like the C u or 
C p is induced by the many-body correlation contributions. Though the mechanism is different 
from the present study, the correlation effect beyond the mean-field approximation may play 
some role in obtaining the coefficients y p and y w . 
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Figure 2: The symmetry energy as a function of density for different models. The MDI(x=0) 
and MDI(x=-l) results are taken from [10]. 

The symmetry energy at high densities from the MDI interaction with x=0 is quite close 
to that given by the SLl*. This is not surprising since the SLl* is rendered to have the same 
symmetry energies as the MDI(x=0) at low densities. On the other hand, this justifies the 
consistency of the symmetry energy given by both models at high densities. One can freely 
adjust y p to fit the symmetry energies given by the MDI(x=-l) model at low densities and then 
examine the behavior at high densities. However, the C p will increase with density by doing 
so, which is contrary to the empirical cases [24, 25]. Therefore, based on the BR scaling the 
symmetry energy at high densities is expected to be between those given by the SLl and the 
SLl*. 

We now turn to the astrophysical implications of the EOS's constrained above. Generally 
speaking, these EOS's have the special features of being soft at moderate densities and stiff at 
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Figure 3: The neutron star mass versus the radius for various models. The parameter set SLl' 
is just the same as the SLl but with y p = 0.526. 

high densities (see Fig. 1). It is thus interesting to compare predictions using these EOS's with 
the recent astrophysical observations. Recent studies on the millisecond pulsar PSR J0751+1807 
suggest that it has a mass of 2.1 ± 0.2(±g:f) with la (2a) confidence [40]. Since the maximum 
mass of neutron star is more sensitive to the EOS at high densities, this requires a rather stiff 
EOS at least at high densities. In Fig. 3, we plot the mass-radius correlation of neutron stars 
obtained from solving the standard TOV equation. In the models with chiral limits, hadron 
masses approach zero at certain critical baryon densities (see Table 1). The maximum central 
energy density in a neutron star is that calculated at the maximum baryon density. With 
SLl, the maximum neutron star mass is 2.04M Q with a radius of R=9.89km. With SL2, these 
two observables are 2.17M0 and R=10.13km, respectively. As the nucleon effective mass at 
normal density is reduced by introducing the parameter y w , the maximum neutron star mass 
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also drops. This is owing to that the introduction of softens the EOS at high densities and 
the energy density at the critical density is thus lowered. With SL1*, the nucleon effective 
mass is M* = 0.6M at po, and the maximum neutron star mass is 1.7M . For a moderate 
value of M* = 0.65M, we can obtain a moderately larger neutron star mass 1.94M . Recent 
measurements on the neutron star EXO 07482-676 gave a mass of M = 2.10 ± 0.28Mq and a 
radius of 13.8 ± 1.8km [41]. The author indicated that the existence of such a massive neutron 
star rules out all the soft EOS's of neutron-star matter. Among the models studied here, 
only the parameter set SL3 can give values close to the measurement for the EXO 07482-676. 
However, the parameter set SL3 gives a pressure at p < 4.6/30 that is too strong compared to 
that constrained by the collective flow data in heavy- ion collisions [6] . For the other parameter 
sets, the radii of maximum mass neutron stars are just around 10km, which is below that given 
in [41], though the maximum mass can be within the error bars of the measurement (except for 
SL1*). 

Compared to the EOS obtained using the FSUGold parameter set [25] the EOS's in the 
present study are stiffer at high densities although they have the same incompressibility k = 230 
MeV at normal density (see Fig.l). The EOS's obtained here thus also result in larger maximum 
masses of neutron stars. In Ref. [5], the authors predicted a radius span of 11.5km<R<13.6km 
for the 1.4M Q canonical neutron star. Using the EOS's obtained in the present work, we obtain 
here a radius span of 11.9km<R<13.1km for the 1.4M Q neutron star, which is comparable with 
that obtained in Ref. [5]. 

4 Summary 

In summary, within the RMF framework we have constructed several new model Lagrangians 
using in-medium hadron properties according to the Brown- Rho scaling due to the chiral symme- 
try restoration at high densities. The model parameters are determined such that the symmetric 
part of the resulting EOS's at supra- normal densities are consistent with that required by the 
collective flow data from high energy heavy-ion reactions, while the resulting density depen- 
dence of the symmetry energy at sub-saturation densities agrees with that extracted from the 
recent isospin diffusion data from intermediate energy heavy- ion reactions. The rearrangement 
terms are found to play an important role in softening the EOS at moderate densities. The 
symmetry energy depends on the in-medium hadron properties that are characteristic of the 
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chiral symmetry restoration at the critical density according to the BR scaling. The resulting 
EOS's are then used to examine global properties of neutron stars. It is found that the radius 
of a 1.4M neutron star is in the range of 11.9km<R<13.1km. Compared to other EOS's, the 
current EOS's have the special feature of being soft at intermediate densities but stiff at high 
densities naturally. This feature is important to produce a heavier maximum neutron star mass 
around 2.0M© consistent with recent observations. 
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